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We present a numerical technique for calculating path integrals in non-compact U(l) and SU(2) 
gauge theories. The gauge fields are represented by a superposition of pseudoparticles of various 
types with their amplitudes and color orientations as degrees of freedom. Applied to Maxwell 
theory this technique results in a potential which is in excellent agreement with the Coulomb po- 
tential. For SU(2) Yang-Mills theory the same technique yields clear evidence of confinement. 
Varying the coupling constant exhibits the same scaling behavior for the string tension, the topo- 
logical susceptibility and the critical temperature while their dimensionless ratios are similar to 
those obtained in lattice calculations. 
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1. Introduction and basic principle 

In this work we present a numerical method for calculating path integrals in Euclidean space- 
time. We apply this method which we call pseudoparticle approach to U(l) and SU(2) gauge 
theories. It is a natural generalization of an idea that was successfully used to show the confining 
properties of ensembles of merons and regular gauge instantons [[I], Qj. 

In the following we explain the basic principle by concentrating on the case of SU(2) Yang- 
Mills theory 1 . We perform Monte-Carlo calculations of SU(2) path integrals 

= - J DA0[A}e- s[A] (1.1) 

F£ v (x) = d^{x)-d v Al{x) + e^ c A^(x)A c v {x) (1.3) 
by representing the gauge field by a linear superposition of a fixed number N of pseudoparticles: 

A a ,(x) = t p ab (iK(x-z(i)). (1.4) 
i=i 

By a pseudoparticle we denote a field configuration p {i)a b n{Xfi —Zp(i)) which has localized ac- 
tion, energy and topological charge (p ab (i): amplitudes and color orientations; Za(i)' positions). 
A common example is the regular gauge instanton, dt{x) = 2rj^ v x v /(x 2 + A 2 ). The integration 
over all gauge fields is replaced by a finite dimensional integral over the amplitudes and color 
orientations of the pseudoparticles and can be computed by Metropolis sampling: 

f N 

DA... - Y\dp(i).... (1.5) 

J i=l 

For certain classes of pseudoparticles it can be shown that in the "continuum limit" (the limit 
of infinitely many pseudoparticles) every gauge field can be expanded as in (L4) [ph. That is in this 



limit the pseudoparticle approach is equivalent to the "full" SU(2) Yang-Mills theory. Therefore we 
expect to get a good approximation of SU(2) Yang-Mills theory even when using a finite number 
of properly chosen pseudoparticles 400 for the results in the following sections). 

In our calculations Euclidean spacetime is restricted to a four-dimensional hypersphere. In- 
side the hypersphere the positions of the pseudoparticles are chosen randomly. Observables are 
measured sufficiently far away from the border. 

2. Maxwell theory 

To check our method we applied it to Maxwell theory. We used Gaussian localized pseudopar- 
ticles 



The corresponding formulas for the U(l) case are very similar. Most of them arise by dropping color indices. 
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and antipseudoparticles 

W = n a ^ve- x2/{nl) (2.2) 

(2-3) 

which form a basis of all transverse gauge fields in the continuum limit: 

A M W = Ep 8 (iK(*-^)) + EW)^(^-^)) ( 2 - 4 ) 

(^v = £ ^v + 8an§ov - $av§on and f}£ v = e a/JV - 5 ajU 5o v + 5 fl vfyu 316 tne t'Hooft- and anti-t'Hooft- 
symbol). There is no need to include longitudinal gauge fields because these fields neither appear in 
the field strength nor in the action. The integration over all gauge fields is replaced by an integration 
over the amplitudes of the pseudoparticles: 

J DA ... = J (rpP^)! flJ^C/M .... (2.5) 

We placed 250 pseudoparticles and 250 antipseudoparticles in a spacetime hypersphere of 
volume 500.0. The average number density per unit volume in this ensemble is n = 1.0 and the 
average distance between neighboring pseudoparticles is d = rT 1 ^ = 1.0. To be able to cover the 
whole spacetime hypersphere with pseudoparticles and to resolve as many ultraviolet details as 
possible the size of the pseudoparticles was set to A = d = 1 .0. 

The potential between two static quarks has been obtained by calculating rectangular Wilson 
loops and using a method [0] which is based on the well known formula 



V q g{R)T + C^-]n{W {RJ) ). (2.6) 

The numerical results and the analytically known Coulomb potential V q q(R) = —e 2 /4nR are plotted 
in Figure [lj For distances R>X = d they are in excellent agreement (A = (/determines the minimal 
size of the field fluctuations; it therefore acts as an ultraviolet cutoff and hence has a similar role as 
the lattice spacing in lattice calculations). 



N - 500, n = 1.0,8= 1.0, \= 1.0 



1 1 



I 



X 



Figure 1: Comparison of the potential of two static charges ±e as a function of their distance obtained 
with the pseudoparticle approach (red points with errorbars) and the Coulomb potential V qq (R) = —e 2 /AnR 
(green curve). 
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3. SU(2) Yang-Mills theory 

For SU(2) Yang-Mills theory we used 1 /r-decreasing pseudoparticles because it has already 
been shown that pseudoparticles with long range interactions are a good choice of parameterizing 
the Yang-Mills path integral ^ 0]. We used three types of pseudoparticles: 

oJW = tfv^xz ("instanton") (3.1) 
tfi( x ) = Wjj.v 2 X +\2 (" ant ii nstanton ") (3-2) 
W = ^^hi (akyron)- (3-3) 



(3.1) and (3^2) are similar to an instanton and an antiinstanton and the akyron 2 (33) is needed to 
get a complete set of functions in the continuum limit. As before the gauge field is a sum over 
pseudoparticles 

A«(*) = £p"^/)a^(x- Z (/))+£p^( 7 >^(^- Z ( 7 0)+Ip flfo W4(x-zW) (3.4) 

i J k 

where p ab {i) is the product of the amplitude a(i) and the color orientation u (i) of the i-th pseu- 
doparticle, i. e. p ab (i) = a(i)u ab (i). The color orientations are restricted by 

u ab (i) = 8 ab (h (i) 2 -h(i) 2 ) +2h a (i)h b (i) + 2e abc h (i)h c (i) (3.5) 



with ho(i) 2 + h(/) 2 = 1. Analogous formulas hold for p ab and p ab . The integration over all gauge 
fields is replaced by an integration over the amplitudes and color orientations of the pseudoparticles: 



J DA... = J (^[\da{i)dh{i)^ \ Y[da(j)dh(j)^ \ Y\da{k)dh{k)^ .... (3.6) 

(dh(i), dh(j) and dh(k) are invariant integration measures on S 3 ). 

A crucial check of the pseudoparticle approach is to verify whether there is confinement or 
not. For that reason we calculated the expectation values of rectangular Wilson loops (W^r)) 
with a fixed ratio of sides R/T = 1/2. We used 400 pseudoparticles (ratio between "instantons", 
"antiinstantons" and akyrons 3:3:2) with a number density per unit volume of n = 1.0 and a size 
of A = 0.5. The results for different coupling constants are plotted in Figure |2|a. The fact that there 
is an area law for large Wilson loops is a clear indication of confinement. We got numerical values 
for the string tension a by fitting an area perimeter function 

- In (w {RJ) \ « a + b(R + T) + oRT (3.7) 

in the usual way. 

The pseudoparticle approach offers the possibility to characterize the importance of certain 
field configurations with regard to confinement. For instance we can turn off the instanton-like 

2 Ancient Greek: akyros = pure gauge (literally "without effect"). 
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N = 400, n = 1 .0, X = 0.5, (# inslanlons) : (# antiinstantons| : {# akyrons) = 3:3:2 N = 400, n = t.0, A= 0.5. g - 4.0 




Figure 2: a) Negative logarithms of expectation values of rectangular Wilson loops — ln(W) as a function 
of their area. Different curves correspond to different coupling constants, b) Comparison of the "akyron 
ensemble" (green curve), the "full ensemble" (blue curve) and the "instanton ensemble" (grey curve). 



field configurations (equations (3J) and (^2)) or the akyrons (equation (33)). Figure |p shows 
that the "full ensemble" and the "instanton ensemble" both exhibit area laws while in the "akyron 
ensemble" the perimeter term is dominating. Performing area perimeter fits yields cv^ron = 0.034, 
CTfuii = 0.318 and a instant on = 0.549, that is a akymn : o M i ■ (Wanton ~ 1 : 9 : 16. This together with the 
fact that a superposition of akyrons has vanishing topological charge density everywhere supports 
the expectation that confinement and topological charge are closely related. 

For a quantitative comparison with lattice results we calculated the topological susceptibility 
X = (Q 2 ) /V ■ Figure |3|a shows the variation of the dimensionful quantities x 1 ^ 4 an d <7 1//2 with 
the coupling constant g. By setting the string tension to the physical value a = 4.2/fm 2 we can 
calculate the size of our spacetime hypersphere. For 1.0 < g < 8.0 the diameter ranges between 
0.51 fm and 2.91 fm. In contrast to dimensionful quantities the dimensionless ratio % [ / 4 /o 1 / 2 is 
nearly independent of the coupling constant (Figure |3|b). Its value 0.295 is in qualitative agreement 



with the lattice result XiJaice/ ] 



V4 / 1/2 



lattice 



0.486 ±0.010 



We also calculated the critical temperature of the confinement deconfinement phase transition. 



N = 400, n - 1 .0, X - 0.5, (# inslanlons) : (# antiinstantons) : (# akyrons) = 3:3:2 



N = 400, o = 1-0,X= 0.5, (# instantons) : (# 



: (# akyrons) -3:3:2 



.*...-....,™*... < „.. < .„ + „^_^I^ 



Figure 3: a) Dimensionful quantities j 1//4 (green curve) and (7 1 / 2 (blue curve) as functions of the coupling 
constant, b) Dimensionless ratio ^ 1 / 4 /<T 1/ ' 2 as a function of the coupling constant. 
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Figure 4: a) Expectation value of the Polyakov loop (L) as a function of the temperature. Different curves 
correspond to different coupling constants, b) Dimensionless ratio T cr - lt i ca \/a 1 ^ 2 as a function of the coupling 
constant. 



To this end we computed the temperature dependence of the expectation value of the Polyakov 
loop which serves as an order parameter. The results for different coupling constants are plotted in 
Figure |]a. Since we consider a finite system and since our ensemble is only approximately center 
symmetric we can not expect to see an exact phase transition. We define the "critical temperature" 
^critical to be that temperature where the expectation value of the Polyakov loop drops below a cer- 
tain threshold, i. e. (L)(r cr j t i C ai) = 1 /2. As shown in Figure |]b the dimensionless ratio ^critical/*? 1 / 2 
is nearly independent of the coupling constant and its value 0.516 is similar to the lattice value 

^critical, lattice/a/^ =0.694 ±0.018 [g. 

4. Conclusions 

We have shown that the pseudoparticle approach applied to Maxwell theory reproduces the 
Coulomb potential for static charges. At the same time the pseudoparticle approach yields an 
area law for Wilson loops in SU(2) Yang-Mills theory. Our results are in qualitative agreement 
with lattice data. We have investigated which properties of field configurations are responsible for 
confinement. Our findings support the common view that topological charge and confinement are 
closely related. 
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